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Abstract 

In the context of random multiplicative energy cascade processes, we derive an- 
' alytical expressions for translationally invariant one- and two-point cumulants in 

logarithmic field amplitudes. Such cumulants make it possible to distinguish be- 
tween hitherto equally successful cascade generator models and hence supplement 
I/-) ' lowest-order multifractal scaling exponents and multiplier distributions. 

o ' 
o 
o . 

Although the underlying hydrodynamic equations are deterministic, the sta- 
tistical description of fully developed turbulence has by now a long tradition 
[1]. Random multiplicative energy cascade models form a particularly simple 
and robust class of such statistical models. While different theoretical models 
can reproduce experimentally observed multifractal scaling exponents rather 
easily [2], observed multiplier distributions [3,4] eliminate many candidate 
cascade generators. Nevertheless, a number of competing generators remain, 
equally successful in reproducing both scaling exponents and multipliers. To 
make further progress in ferreting out the best cascade generator within this 
approximation, new observables are clearly called for. 

While most experiments have concentrated on measuring statistics in the en- 
ergy dissipation density e, we recently found a complete and analytical solution 
working in In e rather than e itself [5] . In this letter, we show that cumulants in 
In e are analytically calculable even when restoring translational invariance to 
the solutions to emulate the spatial homogeneity of experimental turbulence 
statistics. Cumulants turn out to be powerful tools which for third and fourth 
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order differ not only in magnitude but even in sign for the relevant cascade 
generators and can hence be expected to distinguish between different mod- 
els which are otherwise indistinguishable in terms of observed lowest-order 
multifractal scaling exponents and multiplier distributions. 

In the simplest versions of random multiplicative energy cascade models, en- 
ergy flux densities e are generated as follows: in successive steps j — 1, . . . , J, 
the integral scale L is divided into equal intervals of length lj = = L/2- 7 

and dyadic addresses k = (ki---kj) with hi = or 1. At each step j, the 
energy flux e kl ... kj generates fluxes in the two subintervals via 

£k\---kjkj+\ Qki---kjkj+i ^k\---kj i (1) 



where the random variables qi = q^-kjO and qn = q^-kji for the left and right 
subintervals are drawn from a given cascade-generating probability density 
p{Ql, Qr), independently of other branches and generations of the dyadic tree. 
When after J cascade steps the smallest scale 77 = lj = L/2 J is reached, the 
local amplitudes of the flux density field e t t K \ are interpreted as the energy 
dissipation amplitudes at positions 1 < t(n) = (l + Z)/=i ^j^' 1 ^) ^ 2 J in units 
of 77, which are to be compared to experimental time series converted to one- 
dimensional spatial series by the frozen flow hypothesis. 

We have shown previously that, since the product of multiplicative weights 
£ kik 2 -kj — Uj=i Qki-kj becomes additive on taking the logarithm, \n.e klk2 ... kj = 
Z)/=i l n 9fci-fcj, the multivariate cumulant generating function for \ne has the 
analytical solution [5], 



mZ(A ...o,---,Ai...i) = ln/exp A fcl ... fcj ]ne fcl ... fcj 

\ \fei,...,fej=0 

j 1 

,-_i0) ^ki-kj-il) ) 

j = l kl,...,kj-!=0 



(2) 



where the branching cumulant generating function Q (with arguments A^...^ 
J2k j+1 ,...,kj=o ^fci-fcj) i s ^ ne Mellin transform of the cascade generator, 



Q(\l,\r) = ln 



dq L dq R p(q L ,q R )q X L L q R l 



(3) 



which can often be found analytically. A host of analytical predictions for 
statistics in lne for a given cascade generator follow, starting with multivariate 
cumulants obtained directly from In Z through 



C(ki,«2, • • • ,k„) = (hi£ Kl \ne K2 - ■ - \ne Kn ] 



d n \nZ 



(9A K1 (9A K2 



A=0 



(4) 
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Due to the additivity of InZ in (2), these cumulants in Ine become simple 
sums [5] of same-lineage cumulants c n and splitting cumulants c rjS in lng, 



L 



((\nq L Y(\nq R y 



A L =A H =0 

d r+s Q 



c d\ r L d\ s R 



Ar=A R =0 



(5) 
(6) 



where without loss of generality we have assumed Q(Xl, ^r) to be symmetric 
in its arguments. When all n addresses are the same, the n-th order theoretical 
one-point cumulant is simply 



C n («) = ((\ne K ) n ) c = Jc n , 



(7) 



while the theoretical two-point cumulant of order (r, s) for bins K\ — {k\ • • • kjkj + i ■ ■ ■ kj) 
and k 2 = (ki ■ ■ ■ kjk'j +1 ■ ■ ■ k'j) with kj + i ^ kj+i, separated by the ultrametric 
distance D = J — j, is given by [5] 

C r , s {D) = ((\ne K1 Y(\ne K2 y) c = ( J - D)c r+S + (1 - 5 D , )c ryS . (8) 



Three- and higher-point cumulants take on a form very similar to the two- 
point expression [5]. 

Before theoretical cumulants can be compared to experimentally observed 
ones, two complications must be dealt with. The first is that the generat- 
ing function (2) and its cumulants are not translationally invariant, in con- 
flict with the homogeneous statistics characterising experimental results. The 
second complication arises because experimental cumulants are derived from 
measured moments rather than the other way round [6] , requiring translational 
averaging over two-point moments rather than two-point cumulants for theory 
also. The proper procedure is hence to convert theoretical cumulants (8) to 
moments, average these to restore translational invariance, and then convert 
them back to translationally invariant cumulants for experimental comparison. 

For one-point statistics, this procedure is trivial, yielding translationally aver- 
aged one-point cumulants C n = Jc n . Although the latter cannot distinguish 
between a spatially ordered cascade process and one whose field amplitudes are 
spatially randomised, useful information can nevertheless be extracted from 
them. To show this we introduce three factorised model cascade generators, 

P{QL, Qr) = p(Ql) pM , (9) 
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which have proven successful in reproducing multifractal scaling exponents,^ 
namely a simple binomial weight distribution 



Pbinomial(?) 



"'" S{q-{l-a 1 )) + -^S{q-{l + a a )) (10) 



with parameters ot\ = 0.3 and a 2 = 0.65, a lognormal distribution 

-2\ 2 " 



Plognormal (?) 



1 



27rag 



exp 



1 /, 

„ hi q H 

2a 2 1 y 2 



with parameter a = 0.42, and a Beta distribution 

IXA + &) 



PBeta(g) 



r(A)r(/? 2 



(12) 



with parameters j3\ = 4.88 = /?2/7 and q G [0, 8]. All quoted parameter values 
were obtained from fitting to the observed multiplier statistics, including scale 
correlations [4]. 

As shown in Fig. la, all three cascade generators yield almost identical results 
for the lowest-order multifractal scaling exponents r(n) = log 2 (g n ). Since (q) = 
1 by construction, r(l) = for all three distributions. For n = 2 we get 
r(2) = 0.26 for the first two distributions and 0.23 for the beta distribution, 
i.e. indistinguishable within the uncertainty of the experimental intermittency 
exponent fi = 0.25 ± 0.05 [2]. Note that the r(n)'s for the binomial and the 
Beta distributions remain indistinguishable even for n > 3. 

While the multiplier distributions and lowest-order scaling exponents of the 
three cascade generators (10)-(12) are hence demonstrably indistinguishable, 
Fig. lb shows that, while C\ and c 2 are very similar for the three models, the 
same-lineage cumulants exhibit clear-cut differences in third order, with re- 
spective values c 3 = 0.05, 0.00 and —0.05. In fourth order, the same-lineage 
cumulant C4 also exhibits a different sign for the binomial and Beta distribu- 
tion. For the lognormal distribution, C4 is of course again zero. 

The relation between the exponents and cumulants is easily demonstrated 
using the cumulant branching generating function (3): 



d n Q(X,0) 



d n \n(q x ) 



A=0 



In 2 



d n r(\) 



A=0 



(13) 



A=0 



1 We note that other popular log-stable or log-Poisson weight distributions are 
good at reproducing multifractal scaling exponents but fail to account for the proper 
scale-correlations observed in multiplier distributions [4]. 
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i.e. the cumulant c n in In q is the n-th derivative of the scaling exponent r(A) 
at A=0. In principle, therefore, the complete set of same-lineage cumulants 
c n with l<n<oo contains the same information as the complete set of mul- 
tifractal exponents r(n) with 0<n<oo: for the former case an inverse Mellin 
transformation recovers the (factorised) cascade generator, while the latter 
uses the inverse Laplace transform. Once, however, we truncate the two sets 
of observables to the lowest orders, l<n<4, they sample different aspects of 
the complete information. This is the reason why the three cascade generators 
given above can be hard to distinguish by means of lowest-order multifractal 
scaling exponents while still exhibiting significant differences for the third- 
and higher-order same-lineage cumulants. 

As mentioned previously, one-point statistics are not sensitive to the nested 
spatial hierarchy of the cascade. By contrast, eq. (8) shows that two-point 
cumulants depend on the cascade geometry through the ultrametric distance 
D. To access such spatial information, we now concentrate on two-point cu- 
mulants. 

For two-point moments, spatial homogeneity can be emulated by creating a 
theoretical time series consisting of a chain of m — > oo independent cascade 
fields with L/rj = 2 J finest-scale bins each [7]. While this scheme is simple, 
it has previously played a decisive role in explaining scale-correlations for ob- 
served multiplier distributions [3,4] and observed Markov properties [8,9]. For 
two-point statistics with constant bin-bin distance d < L/rj, the appropriate 
averaging is given by 

(m-l)2 J 

PrAd) = (m _ 1)2J £ (^e t y(lne t+d y , (14) 

where t and t+d, corresponding to K\ and k 2 of (8), are in base-ten notation. 
Operationally, these bins are moved over the series in bin-sized steps, suc- 
cessively "seeing" parts of adjacent cascade configurations. For d < L/rj, the 
t-average over the infinitely long chain of the independently and identically 
distributed cascade configurations can hence be replaced by a combination 
of an ensemble average () and a t- average restricted to two adjacent cascade 
configurations, so that (14) simplifies to p rjS (d) = 2~ J J2t=i Pr,s(t, t + d), where 
p rtS (t,t+d) = ((\ne t ) r (\.ne t +d) s )- Since furthermore the two-point correlation 
density factorises whenever bins t and t + d belong to independent cascade 
field configurations, the averaged two-point moment finally becomes 

j (2 J ~d 2 J } 

Pr,s( d ) = E PrAt, t + d)+ E Pr(t)Ps(t + d) . (15) 

z y t=i t=2-'-d+i ) 

Analytic expressions for p r;S (t,t+d) and p r (t) = ((\ne t ) r ) are readily derived 
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by inserting the cumulants (8) into the usual relations between n-variate mo- 
ments and cumulants [6] and thence into (15). 

Because they are functions of the two-point cumulants, two-point densities 
with d > 1 also depend on the ultrametric distance D(t,t + d) between bins t 
and t + d, so that (15) will contain sums of the type 

G n ( J, d) = 1 (J - D(t, t + d)) n . (16) 
z t=i 



In terms of a discrete probability distribution p(D\d, J), proportional to the 
times each value of the D is taken on while t runs its course, this becomes 

G B (J,d)= 5>(£>|d,J)(J-£>r, (17) 

D=l 



with p(D\d, J) empirically found to be 





p(D\d,J) = { 



1 - (d/2 A ) 
d/2 D 



(1 < D < A) 
(D = A) 
(A<D<J) , 



(18) 



where A = |~log 2 d] is the ceiling of log 2 d. Insertion of (18) into (17) leads to 
analytical expressions for the geometrical coefficients 



G (J,d) = (l-2- J d) , (19) 
Gi(J, d) — (J — A) — 2d(2~ A - 2~ J ) , (20) 
G 2 (J, d) = {J - A) 2 - Ad{J - A)2- A + Qd{2- A - 2~ J ) , (21) 

which in turn yield analytical results for the averaged two-point densities 
p r s (d). From these, spatially homogeneous two-point cumulants are constructed 
via the inversion formulae [6] 



C hl (d) 


= Pi,i(d) 


-pi , 




C 2)1 (d) 


= P2,i(d) 


-2piPi,i( d ) 


- P2P\ + 2 pI , 


CzM 


= p3,i(d) 


- 3piP 2 ,i( d ) 


- PsPi - 3p 2 Pi,i( rf ) + Q P 2 iPi,i(d) + 6p 2 p? 


C 2 , 2 (d) 


= p2,2(d) 


-4piP 2 ,i(rf) 


- 2(p lil (d)) 2 -f 2 + Spjp^id) + Ap 2 pl - 



(22) 

3^4 



5^4 



With Eqs. (15)-(22) we arrive for C r ,i(d) at 

C r ,i(d) = Gi(J, d) c r+1 + G (J, d) c.,1 . (23) 
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This turns out to be equivalent to direct translational averaging of (8), i.e. 
2 _J J2t=i Cr,i(t, t + d). For s^l, however, such direct averaging is wrong and 
the full conversion from cumulant to moment to averaged moment and back 
to averaged cumulant is mandatory. For r=s=2 we get, for example, 

C 2 , 2 (d) = d( J, d) (c 4 + 4c 2 c 1;1 ) + G (J, d) (c 2 , 2 + 

-2 [Gi(J, d)c 2 + G (J, d)c 1A } 2 + 2G 2 ( J, d)c\ , (24) 

where the additional terms are a consequence of the quadratic (p x ^cf)) 2 term 
in the expression for C 2)2 (of) in (22). 

Equation (23) has a simple structure, consisting of two terms each of which is 
the product of a geometrical prefactor and a branching cumulant. We note that 
the geometrical prefactor G n ( J, d) depends only on the geometric structure of 
the cascade but not on the cascade generator p(qL, <?r)> in particular, as G n is 
independent of the cumulant order r, it is the same for all C r ^. The cascade 
generator, on the other hand, enters only via the branching cumulants c r+ i 
and c r i. 

For factorised cascade generators (9), eq. (23) simplifies even further to C r . : i(d) = 
G\( J, d) c r+ ±. The distinguishability of c r+ i's for the three model generators 
(10)-(12) therefore implies that two-point cumulants C r ^i(d) of the models 
will differ significantly also as they amplify c r+ i by the geometrical prefactor. 

In order to check whether our analytical results remain statistically significant 
for finite data samples, we simulated a chain of m = 10 7 cascade configurations 
mimicking a time series of m integral lengths, corresponding in size to a typical 
experimental data set. The length of the inertial cascade range was set to 
L/i] = 2 s . One- and two-point densities were sampled according to Eq. (14) 
and then converted into two-point cumulants using relations (22). Variances on 
translationally averaged moments were estimated for each d by picking, with 
a random t, just one pair of energy fluxes (et, £t+d) per two adjacent cascade 
configurations, thereby avoiding correlations between picked pairs. Variances 
for Cr^ were calculated from (22) using standard error propagation. 

Results for C 2 ,i are depicted in Fig. 2, with the solid lines representing the an- 
alytical function Gi(J, d) C3 and simulations yielding the shaded bands. Based 
on Fig. 2, it would appear that statistically significant values for same-lineage 
cumulants c r+ i up to at least r = 2 can likely be extracted experimentally. 
As standard errors increase with increasing order, meaningful statements on 
higher orders become successively more difficult. 

While one-point cumulants C n should no doubt be measured, the two-point 
cumulants contain significant additional information. 
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Firstly, one-point cumulants cannot distinguish between the ordered chain of 
cascade configurations and its spatially randomised equivalent, while two-point 
cumulants can do so. This is true because for the randomised case two-point 
cumulants are simply zero while in the ordered case they retain memory of the 
ordered cascade tree through the geometric coefficients G n ( J, d), which accord- 
ing to relation (17) can be understood as ra-th order moments of the probabil- 
ity distribution p(D\d, J). Thus, two-point cumulants yield information on the 
cascade generator via the branching cumulants c r+ i as well as testing the spa- 
tially nested cascade hierarchy via the geometric coefficients G n (J, d). To date, 
the probing of the treelike structure of the underlying process has not received 
much attention in the literature; we only know of Refs. [10,11], which discuss 
multifractal phase-transition-like behaviour for two-point densities, and of Ref. 
[12], which concerns itself with wavelet scale-scale correlations. 

Secondly, when cascade generators p(qL, qn) do not factorise as in (9), the split- 
ting cumulant c TyS is nonzero, and so same-lineage cumulants plus multifractal 
scaling exponents are insufficient for their complete reconstruction. Indeed, 
more realistic cascade generators, attempting as best possible to reproduce 
the observed one-dimensional energy dissipation data, will typically exhibit 
small but nonzero residual correlations between q L and q^. This expectation 
follows from a projection argument [4] stating that the energy cascade evolves 
in three dimensions but is observed only in one. Through their dependence on 
a nonzero c rjS , the C r>s should thus provide a more complete characterisation 
of the cascade generator. 

To conclude: In the pursuit of finer statistical facets of the energy dissipa- 
tion field in fully developed turbulence, translationally averaged cumulants 
of logarithmic field amplitudes appear to be a promising new tool. We have 
demonstrated for random multiplicative cascade processes that two-point cu- 
mulants can be written as simple products of geometrical coefficients times 
cumulant moments of the cascade generator and can thus distinguish between 
cascade generators which have more or less identical lowest-order multifractal 
scaling exponents and multiplier distributions. Through their dependence on 
the cascade geometry, two-point cumulants are also able to test the treelike 
cascade structure. 

It is tempting to apply two-point cumulants of hie directly to the experimen- 
tal energy dissipation field deduced from hot-wire time series to learn about 
the best approximate random multiplicative cascade model and then to study 
possible dependences on the Reynolds number and the flow configuration. We 
suggest, however, that studies of other and more elaborate models such as 
hierarchical shell models [11], effects of finite inertial range etc. be under- 
taken before this new tool is applied to data. Also, practical problems such 
as choosing the finest resolution scale for the analysis and finite sample size 
effects should be pondered in more detail. 
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Figure 1: (a) Multifractal scaling exponents and (b) same-lineage cumulants 
for the binomial (circles), lognormal (squares), and beta (diamonds) cascade 
generators (10)-(12). 
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Figure 2: Analytical (solid lines) and simulated (shaded bands) two-point cu- 
mulant 6*2,1(0?) as a function of bin-bin distance d, for the binomial (top), 
lognormal (middle) and beta (bottom) cascade generators (10)— (12) . 
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